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We study the interplay between an isostructural critical point and dislocation mediated two- 
dimensional melting, using a combination of Landau and continuum elasticity theory. If dislocations 
are excluded, coupling to the elastic degrees of freedom leads to mean field critical exponents. When 
dislocations are allowed, modified elastic constants lead to a new phase buried in the solid state 
near the critical point. Consistent with a proposal of Bladon and Frenkel, we find an intervening 
hexatic phase near the critical point of the first-order isostructural transition line. Very close to 
the critical point, a transition from the hexatic phase to a modulated bond angle state is possible. 
Ising transitions in symmetrically confined two dimensional colloidal crystals with attractive wall 
potentials require a different model free energy, which we also discuss. 
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' First order phase transitions preserving symmetry and possessing a critical point have been well understood for 
^vq many years, the most common example being the discontinuous liquid-vapor phase transition. Such systems can be 
brought continuously from one phase to another by taking a thermodynamic path which passes around the critical 
^-H , point along the equation of state surface without sudden changes in density. Near the critical point, liquid- vapor 
^ ' systems show critical opalescence among other properties peculiar to second order phase transitions. 

, Recently, attention has been focussed on solid-solid phase transitions . Especially interesting are isostructural 
transitions in which the lattice spacing jumps discontinuously, but the symmetry of the lattice remains unchanged. 
A schematic phase diagram as a function of temperature, T, and particle density, p, is shown in Figure ^. As in 
the liquid- vapor transition, isostructural transitions can have a critical point above which the distinction between the 
l/^ ^ isostructural phases and 5*2 vanishes. Because the solid-solid coexistence region is on the higher density (solid) 
On . side of the fluid-solid coexistence region, a fluid-solid-solid triple point becomes possible. 

Experimentally, bulk isostructural transitions have been observed in high pressure studies of Zr. Above the fee— > uj 
(at 6.7 GPa) and lj — + bcc (at 33GPa) transitions, a lattice symmetry preserving transition between two bee phases 
occurs at 53GPa Q . Another experimentally convenient system with which to study solid-solid transitions is colloidal 
particles. Colloids are often observed to be in an ordered fee or bcc lattice, useful for example, in colloidal array filters 
. The interaction potential between colloid particles can be chemically tailored to control the range of interactions. 
Q For example, one can vary the length of grafted polymers, or change solvent properties Single layers of colloids 
Q : can be confined to study the transitions of 2D solids. This has been achieved by squeezing colloidal particles between 
two flat plates ||^ or floating them at an air- water interface |9| . 
• I— I The theoretical study of the statistical mechanics of colloidal particles has utilized various analytical PJl0|-|l^ and 

numerical methods. It is well established that hard spheres, with or without a superimposed attractive poten- 

. tial, can undergo a solid-fluid phase transition. When the attractive potentials are short-ranged enough, experiments 
- - 1 and theory demonstrate that the liquid-gas coexistence vanishes |l4,n5[. The short-ranged attractive potentials re- 



quired can be achieved with colloidal particles or in computer simulations. A recent numerical study using particles 
with very short-ranged attractive potentials (attractive region < 7% of the hard core diameter), and at high densities, 
show solid-solid phase coexistence. This behavior is consistent with the variational treatment of Tejero et. al. [|j 

In this paper, we explore how ideas of defect mediated melting [ p^ may be incorporated into the schematic phase 
diagram of particles with short-ranged attractions. There is some evidence for a two-stage dislocation mediated 
melting process in liquid crystals colloidal dispersions and magnetic bubble domains |l8| , however, a first 
order transition from solid to liquid (for example by generation of grain boundaries) is also possible JTot . The study 
of defects away from an inherent solid-fluid transition, deep in the solid phase, near a solid-solid critical point, is a 
new opportunity to test predictions of dislocation melting theory [p^ . 

As pointed out by Bladon and Frenkel ||], an intervening hexatic phase (the product of the first stage in the 
two-stage melting scenario) is likely near a two-dimensional isostructural critical point. Chen et. al., studying an 
isothermal, isobaric Lennard-Jones system, point out the large system sizes required in searching for a metastable 
hexatic phase Nevertheless, simulation measurements of elastic constants by Bladon and Frenkel |2| show that 
a dislocation instability (presumably leading to a hexatic phase) must develop in the vicinity of a solid-solid critical 
point. Since the compressibility diverges near the critical point in Figure |l], the effective bulk modulus must soften 
This softening of the bulk modulus promotes dislocation unbinding which signals the solid-hexatic transition. In this 
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paper, we explore these ideas using a simple physical model and determine how the existence of a solid-solid critical 
point may induce a nearby hexatic phase. The effects of elastic degrees of freedom on the Ising critical properties in 
the absence of dislocations are also considered. 

Finally, even though most simulations exhibiting a solid-solid phase transition rely on a very short-ranged attractive 
potential, longer-ranged, soft potentials may also show similar behavior. For example, the three-dimensional solid- 
solid transition in dense Cs and Ce |^ may be pressure induced ||l|,|2^ . Two-dimensional realizations such as confined 
colloids 1^ and colloids or Cgo particles at air-liquid interfaces also possess complicated short and long-ranged 
interactions. However, the two-dimensional melting outlined relies only on a solid-solid critical point and is not 
sensitive to the microscopic mechanism of particle interactions. 

In the next section, we introduce two free energy functionals including both elastic degrees of freedom and the 
Ising-like order parameter describing the isostructural solid-solid phase transition. Elastic deformations and the Ising 
order parameter are coupled via an interaction energy. We also briefly discuss symmetrically confined colloidal crystals 
in 2D with attractive wall potentials. If the wall potential leads to an Ising-like symmetry-breaking, such that the 
crystal prefers to be near one wall or the other, a new type of isostructural critical point becomes possible. 

Section III discusses the simplest model in the absence of defects. Above the critical point, the Ising degrees of 
freedom renormalize the elastic coefficients yielding an effective bulk modulus. Thermal fiuctuations in the (single- 
valued) strains renormalize the critical temperature of the Ising degrees of freedom. We show that the critical 
temperature renormalization is different for the uniform and finite wavevector modes in a manner such that the Ising 
critical point acquires mean-field behavior. The uniform, zero wavevector strains induce an infinite-ranged interaction 
that insures classical critical exponents in this two-dimensional system. 

In Section IV, the strains and Ising order parameter in the presence of a single dislocation are studied by solving 
the extremal (Euler) equations. The strain field configurations are asymptotically identical to those of an uncoupled 
model except they are functions of effective elastic coefficients. The related Ising order parameter near the dislocation 
is also determined. An energy-entropy argument for an isolated dislocation leads to the melting criterion simulated 
by Bladon and Frenkel |^ . Qualitative phase diagrams containing solid and hexatic phases are derived based on this 
criterion. 

The free energy of the full interacting model is expressed in terms of a collection of many dislocations in Section V. 
Renormalization group flows for the the parameters in the presence of a few dislocations are also derived. Dislocation- 
renormalized elastic moduli yield a melting criterion agreeing with simple energy-entropy considerations. The critical 



exponents of the crystal-hexatic transition are the same as for conventional dislocation- mediated melting. |28 



Section VI explores transitions of the hexatic phase induced by the critical point. A modulated hexatic phase is 
possible very near the critical point. If there are no modulated hexatic instabilities, the coupling of the Ising order 
parameter to the hexatic bond angle field leaves the nontrivial Onsager-like critical behavior of the 2D critical point 
intact. 



II. FREE ENERGY MODELS 



In this section we describe the free energy functional for the physics near an isostructural critical point in terms 
of the elastic strain fields, an Ising-like structural order parameter, and an interaction between these degrees of 
freedom. Dislocations, leading to multiple- valued displacement fields will only be considered in detail in later sections. 
Symmetries dictate the form of possible coupling terms. 

A first order transition in a scalar quantity such as density near a critical point can be described by an Ising model. 
In the continuum representation, the "Landau-Ginzburg- Wilson" expression for the free energy takes the form p5[ 

f j2 ( 1 |v7 |2 I 2 I ^ 4 



'^i^ktf^ y (,2'^^' + 2'^ + - + . . .j , (2.1) 

where for a first order liquid-gas transition or isostructural solid-solid transition, if is defined as the density deviation 
from the critical density; ip — p — pc- Near criticality, u is approximately independent of T and r a2{T — Tc)/Tc. 
The linear term proportional to h represents a chemical potential and allows us to study the transitions at densities 
away from the critical isochore. Positive h favors (p > (i.e. the dense phase 52 in Fig. 1) while h < favors iy9 < 



(the dilute solid phase 5*1 in Fig. 1). A potential third order coupling, {wip^, say) in Equation (2.1) can be eliminated 
by a shift in if. 

In addition to an Ising order parameter describing the phase separation, a key feature of the solid phases is their 
resistance to deformations. The elastic deformation energy is p^ ], 

Wei ^ ^ - ^ / d'x (2^4 + Xel) (2.2) 
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where e, 



is the symmetric strain te nsor and /i and A are the Lame coefBcients defined as the micro- 



scopic Lame coefficients divided by fc^T. 
small strains. 



Equation (2.2) is valid for two-dimensional triangular crystals undergoing 



The free energy (2.2) in general contains both smoothly varying single- valued strains as well as singular contributions 
from point-like defects such as dislocations. When a solid is thermally excited, defect co mp lexions must be explicitly 



included in the statistical weight 



as well as in a coupling term, such as Eq. ( p.3| ) below. As discussed in 



Sections IV and V, dislocations can be incorporated by considering the strains as comprised of a smooth "phonon" 
component, Uij, and a singular part, Wy , due to the presence of defects, so that = Uij + Wij. 

We now consider the coupling between the order parameter (p and the strain. Since the model is isotropic, symmetry 
requires that to leading order, only the dilatation, eu, couples to the phase variable. The lowest order interaction 
term for a phase separating system is 



Hi. 



xipeu. 



(2.3) 



We assume g > so that a lattice compression (e« < 0) will bias the system toward (p > 0, i.e. toward the 
denser isostructural solid S2 in Figure [T]. Related couplings between an elastic strain field and an auxiliary field have 
been studied in the context of compressible statistical spin models ^0 31 1, quenched random impurities [ ^^ , and 
vacancy /interstitial defects in crystals. 

IJ31 



For compressible magnets with an Ising-like order parameter, the relevant interaction energy is |30 
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Sx ip^ eu- 



(2.4) 



A model of this kind would be appropriate for colloidal crystals symmetrically confined between glass plates with 
attractive wall potentials of the kind shown in Figure |^. Such a double wall potential can induce (for the appropriate 
pair potentials between the colloidal particles) an Ising-like isostructural critical point inside the two-dimensional 
crystalline phase. Below the Ising phase transition, the crystal prefers to sit closer to one wall or the other. If ip{x) 
represents the deviation of the local particle position from the mid-plane, we can use Eq. (2.1) (with h ~ 0) to 
represent the Ising degrees of freedom and Eq. ( ^.21) to model the elastic ones. The sym met rical confinement (we 
neglect gravity) insures that a coupling of the form ( \2.3i ) is impossible, so we are left with (2.4). 

As discussed, for example, by Bergman and Halperin pC| ], this coupling changes the critical exponents of the Ising- 
like transition in the absence of dislocations provided the specific heat exponent a of the pure Ising system is positive. 
Because a — for the 2D Ising model, there is at most a very slow marginal crossover away from pure Ising-like 
behavior for these symmetrically confined crystals. The coupling g' induces a specific heat-like singularity in the 
inverse elastic bulk modulus near the transition. The renormalized inverse bulk modulus can be defined as 



1 

B 
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T-T, 



(2.5) 



where Uu = / d?xuii{x), A is the system area, and B = fi + X. If the specific heat diverges, B = 0, the energy of one 
isolated dislocation will become small, and we again expect dislocation mediated melting close to the critical point, 
similar to the argument in Ref. |^]. However, this specific heat singularity is much weaker than the singularity which 
induces melting for an asymmetrical model, so a hexatic phase may be difficult to observe experimentally. 

In the remainder of this paper, we consider only the linear coupling, (2.3) where the effects are expected to be more 
dramatic. The total free energy is thus. 



(fx 



A 



d^x 



+ — -hip + 



(2.6) 



Figure ^ shows schematically the coexistence of two different density phases of lattice constant and a_. Under a 
uniform area change, the system responds by undergoing both pure dilatational elastic deformations in each of the 
phases as well as interconversion between local regions of phase, 5*1 or S2. The relative partition between these two 
responses i s go verned by g. In Appendix A, we discuss the stability constraints on the coupling constants determined 
from Eq. (2^) expanded to quadratic order. 
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III. DEFECT-FREE LIMIT 



In this section, we study the interplay between the strain and Ising fields in the absence of defects. The strain 
fields Uij are single-valued and can be thermally averaged out in the quadratic theory. We assume free boundary 
conditions in the {xi,X2) - plane and write the Fourier transformed strain in terms of uniform bulk distortions and 
finite wavevector modes IM 



,(0) 



<f#o 

(0) 



qjU^{q))e" 



(3.1) 



In the decomposition ( p.lD , the uniform symmetric strains ulj must be treated separately because these require three 
degrees of freedom, in contrast to the two independent phonon modes which arise when 7^ 0. A similar treatment 
is required in the study of the zero-wavevector excitation in Bosc-Einstcin condensation. 
Similar to (3.1), we decompose Lp{x) as 



^ ipo + 6(p{x). 



(3.2) 



For the linear coupling model (2.3) considered here, the zero wavevector strain modes couple only to the zero wavevec- 
tor Ising field ipo. After integrating out these uniform bulk strains, we find an effective energy for ipQ. Averaging over 
the fiuctuating crystal displacements, {u{q), q^O), however, gives a different effective energy for the nonzero wavevec- 
tor Ising fields 6tp{x). The final result including the quartic coupling is 









2B J 



2 I " 4 u 



d^x 



ip + B) 



(3.3) 
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where 5lp{x) contains only finite wavevector perturbations of the Ising order parameter about its uniform value (/Sq- 

The q = Q strains {uf^) renormalize the (po modes differently than the finite wavevector phonons {u{q)), renormalize 
the (p{qy^O) modes. The critical temperature of (po is shifted according to T* w Tc{l+g^ /a2B). The mass or quadratic 
coefficient of the ipo field (denoted fp), will, upon decreasing temperature, always vanish before the quadratic coupling 
of the finite wavevector modes, (denoted r). Thus, tq < f, where 
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B 



(3.4) 



and, 



r = r — 



u 2 
2^0- 



(3.5) 



Since the Sip{x) modes remain massive near fp ~ 0, they can be integrated out in perturbation theory to give a 
renormalized effective potential (with no gradient term) for the remaining tpo degree of freedom. It is easily checked 
that the ordering of tpo prevents the subsequent ordering of the Sip{x) modes at a lower temperature within our model. 
We conclude that the critical behavior of a 2D Ising model interacting with strain fields according to Eq. ( ^.3| ) is 
described by mean-field exponents. 

The physical reason for mean-field critical behavior is that the uniform strains induce an infinite-ranged interaction 
between the zero wavevector Ising degrees of freedom. To see this, note that the quadratic contributions to Eq. (3.3) 
may be combined to give 



n 



(2) 
eff 



|V^(f)P+ r 



B 



p^{x) 



A{^i + B)B 



X I d^x' ip{x)ip{x'). 



(3.6) 



The first term is the usual continuum representation of a short-ranged interaction, but the second has infinite range, 
induced by the strain coupling. 
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It is also instructive to consider elastic constants as renormalized by fluctuations of the Ising variable. It is easy to 
show that integrating over the slowly varying Lp degrees of freedom leads to an effective Lame coefficient 



A = A ~ g^/r', (3.7) 

where 

r' = r+'^vl (3.8) 

Consider the behavior of the effective elastic constant A near the (mean field) critical temperature T* discussed 
above. Recall that fQ{T*) ~ 0, and let us assume for simplicity that h ~ Q and T < T* . Since r' = r in this case, we 
have r'{T*) = r{T*) = g^/B ^ g^/{^i + A). It follows that 

A(r*) = -fi (3.9) 

so that the effective bulk modulus 

B(T)=fi + X (3.10) 

vanishes as T — > T*. This vanishing of the effective bulk modulus is expected quite generally, since B is related to the 
Ising fluctuations in our model by an equation similar to (p.5|), namely, 



1 _ 1 

^ ~f f f (3-11) 



The strongly divergent mean fiel d Is ing fluctuations, ^ J (Px J cPx' {(f{x)if{x')) ^ (T — T*) ^, insure that B vanishes. 



The reduction in A implied by (3.7) corresponds to a lattice softening, suggesting that dislocations can more easily 
form. 

Linear stability analysis of the defect-free model is given in Appendix A. The conditions of stability are different 
for zero wavevector and finite wave vector instabilities consistent with the renormalized "masses" for q = and q ^ 
fluctuations, fo and f , respectively. 



IV. ISOLATED DISLOCATION 



In this Section, we determine the behavior of ip and the strains with a single dislocation at the origin. We also 
calculate the free energy of this dislocation and apply an energy-entropy balance to predict when dislocations are 
favored. The dislocation contribution to the field configurations is found by solving the extremal equations developed 
below. 

We first write the total strain as 



- JO) 



+ Wi 



(4.1) 



where 
,(0) 



is the q — part and Wij 



are the singular strains due to the dislocation. After minimizing with respect 
to uf'' and expanding ip{x) about a uniform value ipo + Sip{x), we find that the effective free energy becomes (to 



quadratic order in Sip) 



4! 



4 
^0 



hipo 



Sx 



gSipWi, 



(4.2) 



where fp is given by Eq. ( p.4| ) and r' is given by Eq. (^|^). Near a dislocation ip distorts from ip^. We will assume 
that Lp{x oo) — > ipo, where ipo is one of possibly two uniform minima. This ansatz is appropriate in the presence 
of a single dislocation, since domain walls connecting regions with p — ipo and p ^ —ipo would cost an additional 
gradient energy scali ng a s the system size. 

After minimizing (4^) with respect to ipo, Sip{x), and u, the displacement caused by the dislocation, we obtain 



h+—pQ / (PxSip^ (x) — 0, 



(4.3) 
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and, 



V^5<^ + r'Sip + gwu + O [Stp^) = 0, 



/xV^u + A)V(V • u) + gVS(p = -/xz x bS{x), 



(4.4) 



(4.5) 



where u in the last equation is understood to be the lattice displacement due to the dislocation, Wij — ^{dtUj + djUj). 
These strains which satisfy the extremal equations guarantee vanishing of cross terms between and Wij in ( [t.2| ). 
We the treatment in Ref. |26|1 and write the displacement vector as m = mq + where uq assumes the Burger's vector 
when integrated around a dislocation, i.e., V^uo = —z x bS{x). Note also that V-ito = 0. By taking the divergence of 



(4.5) and substituting in the expression for (V • u) from (4.4), we obtain a closed equation for 5ip{x), 

-2 \ 2^ig 



^i + B 



fi + B 



V • z X bS{x). 



(4.6) 



Note that we have retained the gradient contribution to the energ y in Eq. (4.4). The value of ipQ must be determined 
self-consistently; however, for small Sip{x), the correction to fg in (4^) is small, and (po can be approximated by the 
solution to 



(4.7) 



Having determined the boundary conditions for a single dislocation at the origin, we can solve for a small deviation 
of the Ising order parameter Sip(x). With x ~ {xi, X2}, and b — bei, where e is a unit vector in the xi direction, 



S(p{x) 



^gb 



X2 



TTK'^{fl + B) 



1 - 



2k\x\ 



(4.8) 



where 



(4.9) 



M + S 



A plot of 5ip{x) is shown in Figure |^. The peaks and valleys {x2 > and xi <0 respectively) correspond to a missing 
or extra row of particles. We also calculate the dilatation, Eq. (4.4); 



V • U = Wr 



^ibr' 



X2 



l-OiWe 



(N 



ttk'^{^ + B) |fp 

which is asymptotically inversely proportional to S^p. The total displacements can be found from 



-(/i + \y\/wii — g'S/5if — 2fibS{x)e2- 



(4.10) 



(4.11) 



We note that dwii/dx2 and d6ip/dx2 behave as (5-functions at the origin, and find the total displacements around a 
dislocation, 



Ul 



U2 



b 

2^ 



tan 



Br' — g^ X1X2 



b 

'2^ 



nHfi + B) 



£n\x 



K^H + B) |f|2 

Br' - 5^ xj 



K^^i + B) \x\' 



0(g-'=|5|) 
_0(e-«isi) 



(4.12) 



These displacements are asymptotically identical to those of standard elasticity theory |2g], except with A replaced 
by A. 

Evaluating (4.2) to quadratic order using asymptotic expressions for u and Sip gives the total energy of a single 
dislocation of strength b: 



Kb\ [R 



c(a), 



(4.13) 
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where R is the li near sample dimension and c(a) is function of the dislocation core size a. The most singular part of 
the energy ( 4.13| ) has the same form as the standard result except with effects of the coupling Hint contained in an 



effective Young's modulus 



2fi + X ^ ' 



where the effective Lame coefficient for a dislocation is given by 



A -A- 4, (4.15) 



consistent with Eq. (3.7). As discussed in Section III, this reduction of A causes the effective bulk modulus i? = /i + A 
to vanish near the isostructural critical point. 

The total free energy of a dislocation, Gi, also includes entropy describing configurations associated with the (R/a)'^ 
possible positions the dislocation can occupy. Hence, 

Gi/keT ^Heff -2en{R/a). (4.16) 

The temperature above which the entropy favors dislocation formation is given by Gi/ksT w 0, or, 

K « 167r. (4.17) 

This criterion, modified from the Kosterlitz-Thouless criterion of the decoupled theory {4fiB/{p + B) ^ Wtt) deter- 
mines the phase diagram. The onset of dislocation melting occurs when K > K, where 



K 



4m [Br' ~ g^] 
(u + B)r' - g2 

' ^ (4.18) 
4a^(a* + A) 
2^ + A 



= 167r. 



This condition is equivalent to the usual Kosterlitz-Thouless criterion, with the replacement B ^ B = B — g^/r' . 
Because B{T) as T ^ T* , melting is inevitable near the isostructural critical point. To obtain phase diagrams, we 
assume /i and B to be slowly varying functions of the order parameter tp. If we further assume for simplicity that bare 
elastic constants (unrenormalized by either Ising order parameter fluctuations or dislocations) satisfy /i(pc) ~ A(pc), 
the original dislocation mediated melting theory {g = ), pr edicts dislocation unbinding when ii(pc), HPc ) ^ Stt. 
However, in the modified Kosterlitz-Thouless criterion, ( 1.18| ), the effective parameters may be such that (4.18) is 



satisfied even when /x, A > S tt. In the {ro/u^p — pc) plane, the boundary of a hexatic phase near the isostructural 



critical point as found from (4.18) is shown in Figur e pJ 



We use the solution of the cubic equation (4.7) in ( |4.1^ ) to obtain a quantitative picture in the h — T plane. Phase 
boundaries for /i = A = IOtt in the fo/w and h/u plane are plotted in Figure ||. More precise determination of the 
phase boundary would require a first principles calculation of the unrenormalized elastic moduli p and A. This can be 
obtained for example by density functional theory [ p2[ or primitive cell models [Q. Our simple model however, shows 
salient features of 2D melting induced by an isostructural critical point and agrees qualitatively with the simulated 
phase diagram of Bladon and Frenkel . 



V. EFFECTS OF MANY DISLOCATIONS 



When many dislocations are present, the possibility of more complicated structures arises. In particular, dislocations 
may aggregate to form grain boundaries, or extended domains where </? may take on values different from (poj or for 
fo < 0, both domains of +ipo and —^pq phases may coexist. We will assume these complex aggregates do not form. 
Here, we assume a dilute dislocation gas, and demonstrate that the renormalization group recursion relations for the 
effective parameters are unchanged from those of the uncoupled melting theory. 

The stress tensor associated with Ti. = Hee + + Hint is 

(Jij = (m + A)e«% -I- 2p{eij - -euSij) -f 2g(pSij, (5.1) 
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which can be inverted to give 



1 



OuOij —^pOiJ. 



/i + A ' 

Functional minimization with respect to the phonon displacements yields a divergence-free stress, 







(5.2) 



(5.3) 



to which must be added single- valued fluctuating stresses CTy. Since (5.3) holds away from defects, the extremal stress 
can be expressed in terms of a double curl, Eq] 



(5.4) 



For the present problem, the stress tensor for a collection of dislocations at positions {xa} with Burgers' vectors {6q} 
and "Burger's charge density" b{x) = baS{x — Xa) is given by 



47r 



^im^jndmdn / (f X ' bk{x ')ekl{x - X ') ( 



in 



\X — X 



C 



+ eimejndndm / (f x' G{x - X ')\/'^ ip{x ') + a. 

fji -{- A 



(5.5) 



where V'*G(x) = 6{x) and C is a constant of order unity. The functional form of the free energy H = J (PxaijCij +Hi 
becomes 



H = Hi [^] + Ho [wy ] + H^nt [v, u,,] + Hd [b]HZ' ['P, b] 



where Uij is the smoothly varying part of the strain and 



b{xa) ■ b{xi3)ln 



\Xa-Xi3\\ b{Xa) ■ (fg - Xp)b{xf3) ■ {Xg - Xp) 

\Xa - ^/3p 



knT 



(5.6) 



Y^Wo.)?- (5.7) 



The dislocation contribution to the strain, wa , may be written in terms of Burgers' vectors [2^j2^ 

Ka z-b{xa) 'x {x — Xa) 



Wii{x) 



ATT{fl + A) 



E- 



\X Xq 



(5.8) 



which leads to 



gKo 



47r(/Lt + A) 



Z-b{Xa) X{x — Xa) 



(5.9) 



If the fluctuations in Slp{x) are now averaged out, the effec tive model of interacting dislocations takes the form oiHo, 
except with the replacement Kq K, with K giv en by ( 4.1^ ). Al though we have omitted the details, these results 
follow from straightforward generalizations of ( |4.3| ), (4.4), and ( [4.5[ ). For example, (4.5) is generalized as 

(5.10) 



liV^u+ {fi + A)V(V • u) + gWip — mE ^ ^ baS{x — Xa)- 



The energy (5.6) can be analyzed by rescaling dislocation core sizes and distances in a way similar to that discussed 
in [29| ]. Although recursion relations for g and r as well as the elastic coefficients can be found this way, it is simpler 
to construct recursion relations in terms of the effective coupling K. These recursion relations are pq,E9| 



dyji) 
di 



Btt 



y{l) + 27ry2(£)e^W/i6-j„(^(^)/8^) + Q {v\l)) 



(5.11) 
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di ~ T 



Io{K{£)/87T)~-h{K{e)/87T) 



(5.12) 



where K{£) is defined by (4.18) and y{£) is the dislocation fugacity j/ ~ e ^o/ksT ^ K{£) and y{£) are the usual 



scale-dependent renormalization group coupling constants. The 0{y^) term on the right-hand-side of (5.11) is found 
by considering three dislocation configurations |2§|| . Flows of the parameters, as one considers increasing length scales 
£, which have ?/ — s- indicate states with proliferated dislocations, or a hexatic phase. A m elting criterion identical 
with the modified Kosterlitz-Thouless criterion. Equation ( 4.18 ), arises from analysis of ( ^.11 ). Moreover, we conclude 
that the critical properties of standard dislocation melting theory [ p8| , such as an essential singularity in the specific 
heat, and a divergence of the translational correlation length of the form £.+ {T) ^ exp (const. /\T — T.m\'^'^^^^^), are 
unchanged from the result of the usual dislocation-mediated theory. 



VI. HEXATIC MELTING AND DISCLINATION UNBINDING 



In this Section, we examine the stability of the critical-point-induced hexatic phase. Within the hexatic phase, 
when many dislocations are unbound, the free energy is described by slow spatial variations in the bond angle field 
p7| , p8| . The lowest order symmetry allowed coupling between the bond angle 9 and Ising order parameter has the 
form vV9-Vip, and the total free energy is 



(6.1) 



where is the hexatic stiffness constant; Ka{T) — > oo as the freezing transition is approached from the hexatic phase 
psf . We now consider the behavior of the bond angle and density order parameter in the presence of disclinations. 
The hexatic phase predicted from the standard dislocation melting theory will melt into an isotropic fluid phase upon 
further increase in tem perature. Proliferating disclinations cause this first order transition. 

For our model, ( |6.l| ), we find that the int erac tions between disclinations and ip fields decouple. To see this, we first 
write the extremal equations derived from (6.1); 



KaV'^9 + vV^lp = 



(6.2) 



- V^V3 + + -if^ - h- vV'^e = 0. 
6 

Note that Eq. ( |6.2[ ) allows us to define a conjugate bond angle field 

V 

eijdjO EE diO + -r^dnp. 



(6.3) 



(6.4) 



The density order parameter is assumed to be single-valued; however, the bond angle 9 obeys the noncommutivity 
relation implied by a set of disclinations located at positions {xa}, 



€ijdidj9 = m{x) = X! ^a)- 



(6.5) 



The strength of the a singularity is measured by the charge Sa = ±1. Thus, we find 



di9{x) — eij / d x' ■m{x')djG{x — x') — -—di(p{x) 
J Ka 

Using Eq. (|6.6|), we find that the terms coupling the disclination positions to ip{x) vanish and 



36 



^ s{x)s{x')£n^— 
+v I d^x\I9 ■ Wip 



d^x 



(6.6) 



(6.7) 
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where represents the smoothly varying part of the bond angle field. In addition to a possible disclination unbinding 
transitio n, a finite wavevector instability may develop when higher order terms in are included in (6.7). If we 
expand (6.7) to quadratic order about the minimum of and average over the smooth density fluctuations, we find 



36 



s{x)s{x')h 



\X — X 



Ka 
2 



where we have added higher order terms ii^slV^ 



' and iifc|VV20|2 to Eq. 1^), and defined 

,,2 



Kf 



2^0 



(6.8) 



(6.9) 



Either expression, ( |6.7| ) or (3.8), shows that the interaction energy in (6.1) decouples from the disclination configura- 
tions and only modifies the effective stiffness of the (p or 9q fields. 

Note from Eq. ( |6.9| ) that Kb becomes negative near the isostructural critical point at r — ipo ~ 0. For Kb < 0, a 
finite wavevector instability in bond angle (along with density) occurs at 



± 



\Kb\ 
WSK^ 



1 + J1-3KaKc/Kb, 



(6.10) 



provided 1/q* is larger than the lattice or disclination cutoff size. The thin dotted lines in Figures 5 and 6 show the 
possible modulated hexatic phase. Near the critical point, disclination unbinding and melting directly to an isotropic 
fluid is also possible. 

In principle, the solid could melt directly into a modulated hexatic phase via a first order transition. Curves 
delineating possible hexatic-modulated hexatic and solid-modulated hexatic phase boundaries are shown in Figures |^ 
and ||. 

Provided the gradient coefficient coefficient in (6.7) remains positive the critical point retains non-classical 2D Ising 
critical behavior. Coupling to bond angle degrees of freedom does not induce the long-ranged elastic interactions that 
lead to the mean-field critical behavior discussed in Section III. 
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APPENDIX A: STABILITY ANALYSIS 



Bounds on the elastic parameters in Eq. (5J) follow from a stability analysis in the absence of dislocations. We 
set h = for simplicity. Two cases must be treated separately. First, we analyze stability with respect to the two 
components of the displacements u. Then, stability with respect to the three components of zero wavevector strains 
is checked. 

We first write the Fourier transformed Hamiltonian after decomposing the strains into uniform bulk distortions and 
finite wavevector modes as in Reference [p4| . In the decomposition (3.1), the uniform symmetric strains Uj-^'' must 



be treated separately because these require three degrees of freedom, in contrast to the two independent modes when 
q=/=0. First consider the finite q contribution to the quadratic part of Ti. 



where the phonon dynamical matrix is 



<f#0 



D^j{q) = fiq'^S.j + (fi + X)qiqj. (A2) 
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The eigenvalues of the quadratic form (Al) are given by 



Ao = 



A+ = 



(2a. + A)g + r + g ^ 1 ^^^^^ ^ ^^^^ ^ ^ ^ ^^^^ _ ^^^^ ^ ^^^^^^ ^ ^ ^, 



(A3) 



The onset of instabiUty occurs when an eigenvalue vanishes. Thus, the system is stable against long wavelength 
[q 0) modes when ^ > 0, and r > / (2/i + A). 

For uniform distortions, the quadratic terms in the energy can be written in terms of the three independent 
components of the strain tensor and if 



j2 ' 2 

a X —ix> 
2^ 



where we have defined 



X— ^xx ^yy 
^ — ^xy — ^yx 



The four eigenvalues in this case are given by 

1 



A+ = - 



B + r± V(S-0^ + V 



Al = 2fi, A2 - I 



where B = + X. Stability requires ^ > 0, and 



r > g^ifi + X). 



(A4) 



(A5) 



(A6) 



(A7) 



Thus, as in a pure elastic theory, the uniform strains impose a more stringent criteria on stability than do the finite 
wavevector modes. Therefore, as temperature is lowered (r = a2(T — Tc)), the uniform mode (po becomes unstable 
first. 
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FIG. 1. T — p phase diagram of a system exhibiting a solid-solid phase coexistence. (CP) is a critical point, 5i and 5*2 are 
the low and high density solids respectively. 

FIG. 2. Interaction potential U{S) between particles and flat confining walls. The distance 5 along the ordinate measures the 
distance away from the mid-plane 5 = 0.5 between the walls. When the plates are closely spaced U{S) is given by the dashed 
curve corresponding to an Ising free energy above critical temperature. Increasing the plate gap lowers the temperature of the 
Ising-like free energy (solid curve). 



FIG. 3. Coexistence of two isostructural solid phases. Si and 52 with lattice constants o+ and o-respectively. 



FIG. 4. The solution 5ifi for b = —be\. Note that ip diverges at the dislocation, and a cutoff must be imposed before 5ip gets 
too large; however, unlike the phonon displacement field u, 6ip{x) is single-valued. 

FIG. 5. Phase diagram near the critical point of an isostructural solid-solid transition projected on the (fo /u, if = p — p^) 
plane. For jj, — \ = IOtt, a lobe of hcxatic phase is induced near the critical point. The solid hcxatic-solid boundary corresponds 
to g^/u = 2. The hexatic region becomes smaller for smaller /u. A modulated hexatic phase, evaluated using /Kb = l/457r, 
lies within the thin dotted line. The S\ — S2 coexistence line is approximated by mean-field theory. 

FIG. 6. Phase boundaries projected on the fo/u — h/u plane for /j, = X = lOn and g'^/'u = 2.0. The smaller lobe corresponding 
to a modulated hexatic phase is shown by the thin dotted line. A line of discontinuous transitions is represented by the dark 
line along the x-axis and a renormalized mean-field critical point is at the origin. 
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